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The  parameters  in  Lagrange's  equations  of  motion,  mass,  stiffness,  and 
damping  for  a  mathematical  model  having  fewer  degrees  of  freedom  than 
the  linear  elastic  structure  it  represents  may  be  determined  directly 
from  measured  mobility  data  and  the  approximate  natural  frequency 
associated  with  each  mode.  The  natural  frequencies  are  readily  avail¬ 
able  from  response  plots.  Thus,  using  only  impedance-type  test  data 
without  the  use  of  an  intuitive  mathematical  model,  the  equations  of 
motion  for  the  structure  may  be  obtained  —  a  process  referred  to  as 
system  identification. 


In  conjunction  with  the  determination  of  the  aforementioned  parameters, 
th?  eigenvector  or  mode  shape  and  generalized  mass  corresponding  to 
each  natural  frequency  are  also  calculated. 


A  digital  computer  program  was  generated  to  numerically  test  the  system 
identification  theory.  Computer  experiments  were  conducted  to  test  the 
sensitivity  of  the  theory  to  errors  in  input  data. 
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FOREWORD 


The  work  presented  in  this  report  was  performed  by  Kaman 
Aerospace  Corporation  under  Contract  DAAJ02-70-C-0012 
(Task  1F162204AA4301)  for  the  Eustis  Directorate,  U.  S.  Army 
Air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Virginia.  The  program  was  implemented  under  the 
technical  direction  of  Mr.  Joseph  H.  McGarvey  of  the 
Reliability  and  Maintainability  Division*  and  Mr.  Arthur  J. 
Gustafson  of  the  Structures  Division  .**  The  report  is  pre¬ 
sented  in  four  volumes,  each  describing  a  separate  phase  of 
the  basic  theory  of  structural  dynamic  testing  using  impedance 
techniques . 

^Volume  I  presents  the  results  of  an  analytical  and  numerical 
investigation  of  the  practicality  of  system  identification 
using  fewer  measurement  points  than  there  are  degrees  of 
freedom.  The  parameters  in  Lagrange's  equations  of  motion, 
mass,  stiffness,  and  damping  for  a  mathematical  model  having 
fewer  degrees  of  freedom  than  the  linear  elastic  structure  it 
represents  may  be  determined  directly  from  measured  mobility 
data.  Volume  II  describes  the  method  of  system  identification 
wherein  the  necessary  impedance  data  are  experimentally  deter¬ 
mined  by  applying  a  force  excitation  at  a  single  point  on 
the  structure.  Volume  III  presents  a  method  of  determining 
the  free-body  dynamic  responses  from  data  obtained  on  a  con¬ 
strained  structure.  Volume  IV  describes  a  method  of 
obtaining  the  equations  for  the  combination  of  measured 
mobility  matrices  of  a  helicopter  and  its  subsystems.  The 
response  of  the  combination  of  a  helicopter  and  its  sub¬ 
systems  is  determined  from  data  based  on  the  experimental 
results  of  the  main  system  and  subsystems  separately. 


^Division  name  changed  to  Military  Operations  Technology  Division. 
**Division  name  changed  to  Technology  Applications  Division. 


v 


TABLE  OF  CONTENTS 


Page 

ABSTRACT .  iii 

FOREWORD .  v 

LIST  OF  ILLUSTRATIONS .  ix 

LIST  OF  TABLES .  x 

LIST  OF  SYMBOLS .  xi 

INTRODUCTION .  2. 

THEORY .  2 

DERIVATION  .  2 

DERIVATION  OF  THE  DOMINANT  MODE  EIGENVALUE 

PROBLEM .  5 

IDENTIFICATION  OF  STRUCTURAL  PARAMETERS .  9 

PARAMETERS  OF  THE  MATHEMATICAL  MODEL .  11 

ITERATION  PROCEDURE .  12 

INTERPRETATION  OF  ELEMENTS  IN  THE 

REDUCED  MASS  MATRIX .  14 

THE  REDUCED  MASS  MATRIX .  16 

INFORMATION  LOSS  IN  MATRIX  INVERSION .  20 

HOW  TO  MINIMIZE  INFORMATION  LOSS .  23 

WHEN  A  CALCULATION  FAILS .  28 

THE  T  MATRIX  AND  MODAL  PARAMETERS .  29 

IDENTIFIED  GENERALIZED  MASSES .  32 

RESPONSE  FROM  IDENTIFIED  MODEL .  39 

CONCLUSIONS .  46 

LITERATURE  CITED .  47 

vii 


TABLE  OF  CONTENTS  (Continued) 


Page 

APPENDIX  -  COMPUTER  PROGRAM  DESCRIPTION . 48 

DISTRIBUTION  .  138 


viii 


LIST  OF  ILLUSTRATIONS 


Figure  Page 


to  " 

2  Five-Point  Model  Response  Obtained  From 
Equations  With  Identified  Parameters; 

Driving  Point  at  Hub .  40 

3  Nine-Point  Model  Response  Obtained  From 
Equations  With  Identified  Parameters; 

Driving  Point  at  Hub  ......  .  41 

4  Twelve-Point  Model  Response  Obtained  From 
Equations  With  Identified  Parameters; 

Driving  Point  at  Hub .  42 

5  Five-Point  Model  Response,  Effect  of 

Model;  Driving  Point  at  Station  1 .  43 

6  Nine-Point  Model  Response  Effect  of 

Model;  Driving  Point  at  Station  1 .  44 

7  Twelve-Point  Model  Response,  Effect  of 

Model;  Driving  Point  at  Station  3 .  45 


LIST  OF  TABLES 


Table  Page 

I  IDENTIFICATION  OF  GENERALIZED  MASSES, 

5X5  MODEL  OF  20  X  20  SPECIMEN . .  33 

II  IDENTIFICATION  OF  GENERALIZED  MASSES, 

5X5  MODEL  OF  20  X  20  SPECIMEN . .  .  34 

III  IDENTIFICATION  OF  GENERALIZED  MASSES, 

9X9  MODEL  OF  20  X  20  SPECIMEN .  35 

IV  IDENTIFICATION  OF  GENERALIZED  MASSES, 

9X9  MODEL  OF  20  X  20  SPECIMEN .  36 

V  IDENTIFICATION  OF  GENERALIZED  MASSES, 

12  X  12  MODEL  OF  20  X  20  SPECIMEN .  37 

VI  IDENTIFICATION  OF  GENERALIZED  MASSES, 

12  X  12  MODEL  OF  20  X  20  SPECIMEN .  38 


x 


LIST  OF  SYMBOLS 


the  damping  matrix 


a  damping  matrix;  [d]  =  u>[c];  for  damping  forces 
which  are  proportional  to  displacement 


vector  of  external  forces  acting  along  the 
generali2ed  coordinates 

force  phasor,  {f}  =  {f}e^'U)t 


the  structural  damping  coefficient  of  the  i-th  mode 


i  or  j  indices;  imaginary  operator  (i  =  /^T) 


the  generalized  stiffness  of  the  i-th  mode 


the  stiffness  matrix 


the  generalized  mass  of  the  i-th  mode 


the  mass  matrix 


N  or  n  the  number  of  degrees  of  freedom  in  the  structure 


vector  of  velocities  of  the  generalized  coordinates 


velocity  phasor,  {y :  =  {yje"1 


[Y .  .  ]  matrix  of  mobilities  at  forcing  frequency  w; 


IY(U)]=  I9yi/9£j]  (u) 


xi 


LIST  OF  SYMBOLS  (Continued) 


* 

Yi(u))  generalized  mobility  of  the  i-th  mode  at  forcing 
'  '  frequency  to 

[Y]  natrix  of  acceleration  mobilities 

[z^j]  matrix  of  impedances  at  forcing  frequency  to; 

z* .  v  generalized  impedance  of  the  i-th  mode  at  forcing 
1  frequency  to 

_* 

z. .  .  complex  conjugate  of  the  i-th  mode  generalized 
1  impedance  at  forcing  frequency  to 

it 

|z..  v|  absolute  value  of  the  i-th  mode  generalized  impedance 
1 at  forcing  frequency  to 


{y}.  the  i-th  column  of  [T]  ;  the  gamma  vector  of  the  i-th 
1  mode;  a  left-hand  eigenvector  of  [k]-1[m] 

—1  -T 

[ T )  the  left-hand  eigenvectors  of  [  k]  [m] ;  [$] 


6^-*  Kronecker's  delta 

{ 4> } ^  the  modal  vector  of  the  L-th  mode 
[$]  matrix  of  modal  vectors 


to  forcing  frequency 


the  natural  frequency  of  the  i-th  mode 


LIST  OF  SYMBOLS  (Continued) 


SUPERSCRIPTS 

R  the  real  part  of  a  complex  quantity 

I  the  imaginary  part  of  a  complex  quantity 

*  a  generalized  parameter  associated  with  a 

particular  mode 

T  the  transpose 

-T  the  inverse  transpose 

SUBSCRIPTS 

(cd)  the  forcing  frequency  at  which  the  quantity  was 

measured  or  calculated 

k  forcing  frequency 


A  dot  over  a  quantity  indicates  differentiation 
with  respect  to  time 

BRACKETS 

[  1 , (  )  matrix 

["  J  diagonal  matrix 

{  }  column  or  row  vector 


mf*9' 


’WWIWW' 


INTRODUCTION 


The  success  of  a  helicopter  structural  design  is  highly 
dependent  on  the  ability  to  predict  and  control  the  dynamic 
response  of  the  fuselage  and  mechanical  components.  Conven¬ 
tionally,  this  involves  the  formulation  of  intuitively  based 
equations  of  motion.  Ideally,  this  process  would  reduce  the 
physical  structure  to  an  analytical  mathematical  model  which 
would  predict  accurately  the  dynamic  response  characteristics 
of  the  actucl  structure.  Obviously,  the  creation  of  such  an 
intuitive  abstraction  of  a  complicated  real  structure  requires 
considerable  expertise  and  inherently  includes  a  high  degree 
of  uncertainty.  Structural  dynamic  testing  is  required  to 
substantiate  the  analytical  results.  The  analysis  is  modified 
until  successful  correlation  is  obtained  between  the  analytical 
predictions  and  the  test  results.  Finally,  the  mathematical 
model  can  be  used  to  incorporate  changes  to  improve  the  struc¬ 
tural  integrity  of  the  helicopter. 

This  report  describes  the  theory  of  structural  dynamic 
testing  using  impedance  techniques  as  applied  to  a  mathe¬ 
matical  model  having  fewer  degrees  of  freedom  than  the  struc¬ 
ture.  Reference  1  describes  the  method  of  obtaining  a  model 
directly  from  test  measurements  for  a  hypothetical  structure 
which  has  the  same  number  of  degrees  of  freedom  as  the 
mathematical  model.  In  reality,  the  number  of  degrees  of 
freedom  of  a  physical  structure  is  infinite;  therefore,  the 
usefulness  of  model  identification,  necessarily  with  a  finite 
number  of  degrees  of  freedom,  using  impedance  testing  tech¬ 
niques  depends  on  the  ability  to  simulate  the  real  structure 
with  a  small  mathematical  model. 

The  process  of  deriving  the  equations  of  motion  from  test 
data  is  referred  to  as  system  identification.  The  only  input 
information  required  in  this  theory  is  measured  mobility  data 
and  the  approximate  natural  frequency  of  each  mode.  This 
information  can  be  obtained  from  impedance  testing  of  the 
actual  structure  over  the  frequency  range  of  interest  yielding 
the  second  order,  structural]  ly  damped  linear  equations  of 
motion . 

System  identification  theories  to  be  of  any  practical  engi¬ 
neering  significance  must  be  functional  with  a  reasonable 
degree  of  experimental  error.  In  this  report,  a  series  of 
computer  experiments  incorporating  experimental  errors  was 
documented.  This  report  presents  a  modification  and  exten¬ 
sion  of  the  analysis  derived  in  Reference  1  such  that  an 
identified  model  with  a  finite  number  of  degrees  of  freedom 
simulates  the  actual  structure  wherein  the  number  of  degrees 
of  freedom  is  infinite. 


1 


THEORY 


DERIVATION 

The  equations  of  motion  in  matrix  form  of  a  linear  system 
are,  as  shown  in  Reference  1, 


[m]  {y}  +  [c]  {y }  +  [k]{y}  =  { f } 


(1) 


Assume  a  steady-state  solution  of  the  form 

{y}  =  {y}eiwt  and  {f }  =  {f}eiwt 
Substitute  these  equations  into  Equation  (1)  to  give 


l?[m]w  -  £[k])i  +  [cl]  (y>  =  ih  (2) 

or 

(i[z^]  +  [z^]){y}  =  [zj{y}  -  { f } 

where  zij (w)  is  defined  herein  as  the  element  velocity  impedance 
measured  at  w. 


The  element  impedance  can  also  be  expressed  as 

z .  .  .  x  =  9f . /3y . 

13  v<i')  1'  J  j 


— T  T 

If  Equation  (2)  is  premultiplied  by  [$]  [4>]  and  post- 

multiplied  by  [4)]  [4>]  1  where  [4>]  is  the  matrix  of  modal 
vectors,  the  result  is 


m"T[i(mT[m]  [«]u  -  J[4»]T [k]  [#]  )  +  [4>]T[c]  [*]]  [*] 

The  diagonal  generalized  mass  is  expressed  by 

tTJ  =  [$]T[m]  [*] 

The  diagonal  generalized  stiffness  is  given  by 

m  =  mTlk]  [*] 


(3) 

(4) 

(5) 
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Assume  that 


^PgHJ  =  mT[c][*]  (6) 

such  as  would  be  expected  from  structural  damping  in  a  lightly 
damped  structure.  Substituting  Equations  (4),  (5)  and  (6) 
into  Equation  (3)  yields 


t4  (<u)  J 


Define  the  i-th  modal  impedance  as 


•*  .  Xiv  ,  ^i^l 

zi(u)  =  >  +  — 


(7) 


and  substitute  into  Equation  (7)  to  give 

,-Tf 


"  141 


[f-*  lr.,-1 


(8) 


The  elemental  mobility  at  forcing  frequency  w  is  defined  as 

Y.  .  ,  .  5  8y./3f j  and  is  equal  to  the  ratio  of  the  velocity 
13(10)  Jr  j  M 


phasor  along  the  coordinate  i  to  the  external  force  phasor 
along  the  coordinate  j  when  no  other  forces  are  externally 
applied.  The  full  mobility  matrix  is  given  by 


^(003=  t9y/«)(u)-  t3«/3y  1  (w)_1  =  b,.,!’1 


(9) 


Therefore,  using  Equation  (8)  it  is  seen  that 


‘Y(-i>  *  «*« 


1 

(<o)\ 


i*iT  5 


The  modal  mobility  of  the  i-th  mode  measured  at  to  is 


. *  . *p  ,  . *T 

Y.  ,  .  =  Y . ,  .  +  iY.  ,  . 
i(w)  i(to)  l((o) 


*R  .  *1 

Z .  -  iZ.,  . 
1(0))  i(to) 


i  (to) 

to 


.  * 

Zi  (to) 

•  *  2 


. .  _  to- 
-  i  r?iU  ~) 


..*R  *  ,-*I  2  ^i*!,2  +  m  M  -  *L)2 

Ui(U))  +  (Zi(<o))  +  to 


(10) 
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Dividing  numerator  and  denominator  of  the  previous  equation 
by  the  generalized  mass  ^  ^ 


.  * 

yi(o» 


9j*i 
■ni“ 


—  i  (cu  — 


Xi 


.05' 


_  r-  n±- 

/  \  2  .  ^  /  %i  x  2 

^i  ii)^  ^i  w 


Substituting  the  natural  frequency  of  the  i-th  mode 


.  * 

Yi  (u) 


g .  ft .  2  ft .  2 

'W7.  (  -1-  ■■---)  +>>7.  (u>  -  — ) 

05  *1  0) 


Separating  this  equation  into  the  real  and  imaginary 
components  yields 


.  * 

YiU> 


— (— )2 

<V  Y 


-  l 


"I) 

jr 


05 


+  (^y  -1) 


(11) 


Finally,  from  Equation  (iO)  ,  the  real  mobility  may  be 
written  as 


(i2» 

Reference  1  indicated  that  because  the  real  modal  mobilities 
of  modes  far  removed  from  the  forcing  frequency  become 
negligible  compared  to  adjacent  modes,  the  real  mobility 
matrix  at  any  frequency  is  ordinarily  affected  only  by  vnodes 
in  close  proximity  to  the  forcing  frequency.  The  measured 
real  mobility  matrix  at  a  particular  frequency  reflects  the 
influence  of  only  the  most  dominant  modes  in  that  frequency 
of  measurement  region.  Therefore,  it  is  unrealistic  to  use 
the  real  mobility  matrix  measured  at  any  specific  frequency 
to  determine  parameters  other  than  those  associated  with 
neighboring  modes. 
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Reference  1  also  shows  that  the  imaginary  modal  mobilities 
of  modes  associated  with  frequencies  less  than  the  forcing 
frequency  asymptotically  approach  a  constant.  An  imaginary 
mobility  matrix  contains  the  effect  of  all  lower  modes  in 
proportion  to,  or  greater  than,  the  magnitudes  of  their 
generalized  masses.  Therefore,  it  is  impractical  to  use 
imaginary  mobility  matrices  to  evaluate  properties  associated 
with  natural  frequencies  far  above  the  forcing  frequency. 

These  characteristics  of  the  modal  mobility  make  it  impossible 
to  determine  the  system  parameters  from  the  n  equations  in  n 
unknowns  obtained  from  mobility  matrices  measured  at  any  two 
forcing  frequencies . 

Even  if  the  modal  mobility  were  amenable  to  determination  of 
the  system  parameters,  the  precision  of  measurement  which 
would  be  required  to  do  this  for  most  systems  is  impossible 
to  achieve.  The  modal  approach  derived  below  avoids  this 
problem. 


DERIVATION  OF  THE  DOMINANT  MODE  EIGENVALUE  PROBLEM 


Equation  (10)  may  be  written 
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N 

E  Y, 
i=l 
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where  {d>}  is  a  column  in  [$]  and  N  is  the  order  of  the 
matrices.  Define  [T]  =  [$]-T,  and  Equation  (8)  may  be  written 
as 
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where  {y}  is  a  column  in  [T] 
Each  matrix 


ix  Yi(u)^i^^i  an1  ^ — {yHyi)T  in  Equations  (13) 
*“  J  ^-Yi(u)) 


and  (14)  is  of  rank  one,  but  the  summation  of  as  many  of  these 
successive  modal  matrices  as  the  order  N  of  the  matrix  is  a 
nonsingular  matrix. 
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Similarly, 
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The  iteration  procedure  used  to  solve  the  eigenvalue  problem 
in  Reference  1  employed  the  imaginary  part  of  a  mobility 
matrix  measured  at  a  frequency  just  above  the  N-th  natural 
frequency.  The  method  used  to  solve  the  eigenvalue  problem 
in  the  present  report,  which  was  found  to  give  more  accurate 

results,  utilizes  the  sum  of  the  real  parts  of  the  mobility 
matrices  measured  near  each  of  the  natural  frequencies 
associated  with  the  actual  model.  It  has  been  indicated 
previously  that  a  measured  real  mobility  matrix  reflects  the 
influence  of  only  the  most  dominant  modes  in  the  vicinity  of 
the  forcing  frequency.  Therefore,  summation  of  a  discrete 
set  of  the  real  mobility  matrices  measured  at  forcing  fre¬ 
quencies  near  the  corresponding  natural  frequencies  should 
contain  precisely  the  information  relevant  to  the  model 
normal  modes. 


The  eigenvalue  problem  may  be  formulated  as  follows.  Con¬ 
sider  the  summation  of  the  real  mobility  matrices  measured 
at  a  discrete  set  of  frequencies  near  the  first  NR  natural 
frequencies.  Take  the  inverse  of  this  matrix  and  pre¬ 
multiply  by  a  real  mobility  matrix  measured  at  any  frequency 
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If  Equation  (16)  is  postmultiplied  by  there  results 
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but  yields  a  column  matrix  comprised  of  zeroes 

except  for  a  1  in  the  i-th  position.  Finally, 
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The  eigenvalue  problem  is  finally  formulated  as 
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(17) 


If  the  order  of  multiplication  is  reversed  in  Equation  (16)  , 
an  eigenvalue  problem  is  developed  in  which  the  eigenvector 
is  the  gamma  vector  of  the  i-th  mode.  Consider  the  same 
parameters  as  in  Equation  (16);  only  the  order  of  multiplica¬ 
tion  of  the  matrices  is  changed. 
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By  definition,  [ r ]  =  [<t>]  T  and  [i‘]  ^  =  [$]T;  substituting  into 
Equation  (18)  yields 
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If  Equation  (19)  is  postmultiplied  by  iy}-,  a  column  of  [T], 
and  the  same  procedure  is  followed  as  was  used  in  obtaining 
Equation  (17),  Equation  (19)  becomes 
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which  is  an  eigenvalue  problem  with  the  eigenvector  equal  to 
the  gamma  vector  of  the  i-th  mode. 


IDENTIFICATION  OF  STRUCTURAL  PARAMETERS 


A  successful  identification  procedure,  using  normal  mode 
techniques,  should  separate  the  effect  of  each  mode  in  a 
mathematical  sense,  regardless  of  the  number  of  stations 
where  mobility  measurements  are  taken  on  the  structure.  If 
satisfactory  normal  mode  separation  required  a  certain  minimum 
number  of  measurement  stations  greater  than  the  number  of 
degrees  of  freedom  chosen  for  the  model,  the  most  that  can 
be  expected  is  an  approximate  model,  possibly  including 
optimization  procedures  designed  to  satisfy  all  system  con¬ 
straints.  This  situation  is  considered  in  detail  in  Reference 
2  in  which  a  mathematical  model  is  derived  from  test  data 
such  that  identification  of  the  structure  is  obtained  closest 
to  any  specified  analytical  approximation. 

Satisfactory  normal  mode  separation  requires  that  the  values 
of  * *R  and  **I  be  independent  of  the  number  of  degrees  of 
Vj)  Z(u)j) 

freedom  of  the  model.  The  values  of  the  generalized  mass  (^)  , 
the  corresponding  identified  natural  frequency  (^)  ,  and  the 
generalized  stiffness  as  defined  below  are  then  also  in¬ 
dependent  of  the  number  of  measurement  stations. 
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The  two  forcing  frequencies  (cok)  and  (ok)  are  chosen  in  the 

vicinity  of  the  corresponding  natural  frequency  which  is 
available  from  test  data.  The  generalized  impedance  of  the 
i-th  mode  at  forcing  frequency  (w)  is  obtained  from  the 
generalized  mobility  of  the  i-th  mode  at  forcing  frequency 
(oj)  .  It  follows  from  Equation  (13)  that  the  modal  mobilities 
are  given  by 


N*  1 
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and,  therefore,  the  orthogonality  condition  for  gamma  vectors 
is 


The  modal  impedance  of  the  i-th  mode  at  uij  is 
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It  follows  that 
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The  damping  coefficient  for  the  i-th  mode  is  most  readily 
given  by 

*R 

"j2iua> 

9i  *  if.  9  <25> 


which  follows  directly  from  Equation  (7)  .  The  damping  co¬ 
efficient  for  the  i-th  mode  may  also  be  obtained  by 


Using  a  measurement  of 
resonance ,  the  damping 
Equation  (11)  as 
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real  Mobility  taken  precisely 
coefficient  may  be  calculated 


(26) 


at 

using 
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PARAMETERS  OF  THE  MATHEMATICAL  MODEL 


The  elements  of  the  influence  coefficient  matrix,  being  a 
measure  of  displacement  per  unit  force,  are  independent  of 
the  number  of  measurement  stations  defining  the  order  of  the 
matrix.  Conversely,  the  elements  of  the  stiffness  and  mass 
matrices  assume  different  values  as  the  number  of  degrees  of 
freedom  of  the  model  is  changed.  The  identification  procedure 
used  in  Reference  1  calculates  both  stiffness  and  mass 
matrices  by  summing  the  effects  of  each  consecutive  mode 
and  defining  the  incomplete  matrices  as  the  sum  up  to  and 
including  a  particular  mode.  If  the  order  of  the  model 
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degrees  of  freedom  is  changed  from  ND  for  the  structure  to  NR 
for  the  model,  the  corresponding  incomplete  mass  and  stiff¬ 
ness  matrices  will  not  be  directly  comparable,  on  a  modal 
basis,  to  the  structure  mass  and  stiffness  matrices.  It  is 
more  expedient  to  identify  the  influence  coefficient  matrix 
[c]  and  the  inverse  of  the  mass  matrix  [jq] .  Premultiplying 
Equation  (4)  by  [4>]-T  and  postmultiplying  [$]“1  and  taking 
the  inverse  of  the  resulting  equation  yields 

tm  =  [m]  1  =  E  —  (<j>.  }{<}>.}  (28) 

i=l^i  1  1 

If  the  same  operations  are  performed  on  Equation  (5) ,  the 
result  is 

nr 

[c]  =  [K]  1  =  E  {<j>.  }T  (29) 

i=l  Q' lnji 

Set  [c]  =  ~  td ]  and  using  Equation  (6)  there  results 

irsdfi  =  i[tiTtd)[4] 

Solving  for  the  damping  matrix  yields 

[d]  =  4>-T  TgXJ  $-1 

Substituting  [T]  =  [<j>]  T,  [r]T  =  [0]  ^  and  HfJ  =  into 

the  previous  equation  gives 

[d]  =  [T]  ir]T 

The  damping  matrix  can  be  expressed  as 

NR  o  m 

[d]  =  gi^i2/^i{YiHyi}T  (30) 


ITERATION  PROCEDURE 


The  calculation  of  the  modal  parameters  such  as  generalized 
mass,  stiffness  and  the  corresponding  natural  frequency  re¬ 
quires  the  generalized  impedance  at  a  particular  frequency 
for  each  mode  under  consideration.  The  modal  impedance  is 
a  function  of  the  generalized  mobility  for  the  same  mode  and 
forcing  frequency.  As  indicated  in  Equation  (24),  the  modal 
mobilities  are  dependent  upon  the  matrix  of  gamma  vectors  and 
its  transpose.  The  iteration  process  as  originally  formulated 
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in  the  present  work  involved  iteration  on  the  normal  mode 
vectors  with  a  subsequent  inversion  operation  to  determine 
the  gamma  vectors.  This  sequence  introduced  errors  into  the 
system,  with  the  result  that  the  gamma  vectors  did  not  re¬ 
semble  the  assoc  .ated  gamma  vectors  obtained  from  the  specimen 
representing  the  actual  structure.  The  iterated  normal  mode 
vectors  obtained  from  the  mathematical  model  were  extremely 
close,  particularly  at  the  lower  modes,  to  the  specimen,  or 
exact,  normal  mode  vectors.  Nevertheless,  any  discrepancy 
between  the  model  iterated  modal  vectors  and  the  exact  values, 
however  small,  was  magnified  in  the  inversion  process,  causing 
the  gamma  vectors  to  bear  little  resemblance  to  the  specimen 
gamma  vectors.  Therefore,  it  was  deemed  advisable  to  iterate 
on  the  gamma  vectors  directly  and  disperse  with  the  inter¬ 
mediate  inversion  operation. 

To  equalize  the  effect  of  each  modal  mobility  in  the  matrix 
iteration  Equation  (20),  several  normalization  procedures 
were  incorporated  into  the  method.  First,  each  real  mobility 
matrix  was  normalized  on  the  largest  element  of  the  respec¬ 
tive  matrix.  This  procedure  proved  satisfactory  except  in 
some  situations  where  the  elements  of  the  mobility  matrices 
were  approximately  of  the  same  magnitude  but  the  largest 
elements  differed  in  algebraic  sign.  This  resulted  in  a 
cancellation  effect  among  the  real  mobility  matrices  and 
an  incorrect  summation,  thereby  causing  erroneous  calculated 
gamma  vectors.  A  modification  to  the  normalization 
technique  was  applied  whereby  the  real  mobility  matrices  at 
each  forcing  frequency  were  divided  by  the  absolute  value  of 
the  largest  element  in  the  respective  matrices.  As  a  further 
refinement  on  the  normalization  procedure,  the  real  mobility 
matrix  at  each  forcing  frequency  was  normalized  on  the  r  ot 
mean  square  associated  with  each  respective  matrix.  Occasion¬ 
ally,  these  operations  also  caused  problems  in  the  final  modal 
generalized  mass  and  natural  frequency  calculations.  For 
example,  if  a  mobility  matrix  calculated  at  a  particular  fre¬ 
quency  contained  one  element  that  dominated  the  matrix, 
normalization  of  the  mobility  matrix  on  this  element  would 
effectively  submerge  the  influence  of  the  matrix  in  the 
summation  of  the  real  modal  mobilities.  Again,  the  calculated 
modal  parameters  would  obviously  be  incorrect.  Similarly,  if 
several  elements  of  the  mobility  matrix  measured  at  a  specific 
forcing  frequency  were  of  greater  magnitude  than  the  remaining 
elements,  the  root  mean  square  value  would  be  affected  and 
normalization  by  this  value  would  yield  a  matrix  wherein  the 
elements  were  substantially  reduced.  Therefore,  any  such 
matrix  would  not  be  realistically  represented  in  the  summation 
f  the  real  mobility  matrices;  consequently,  the  modal  gen¬ 
eralized  mass  and  natural  frequency  would  be  incorrect. 
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Finally,  each  mobility  matrix  was  multiplied  by  the  respective 
forcing  frequency  yielding  an  acceleration  mobility.  These 
acceleration  mobility  matrices  were  substituted  for  the 
velocity  mobilities  appearing  in  Equation  (17)  and  Equation 
(20)  when  iterating  for  the  mode  shapes  and  gamma  vectors 
respectively.  This  technique  was  also  plagued  with  similar 
difficulties  that  the  aforementioned  normalization  procedures 
incurred.  Fortunately,  when  the  computer  experiments  were 
executed  incorporating  any  of  the  previously  discussed 
normalization  methods,  the  conditions  which  yielded  erroneous 
results  were  readily  discernible.  In  these  instances,  the 
calculations  for  the  modal  generalized  mass  and  natural  fre¬ 
quency  produced  results  which  were  obviously  incorrect. 

For  the  conditions  which  were  recognized  to  be  in  error,  the 
computer  experiments  were  reevaluated  substituting  a  different 
normalization  option.  Generally,  the  results  obtained  by 
altering  a  normalization  procedure  yielded  modal  parameters 
which  were  correct. 

INTERPRETATION  OF  ELEMENTS  IN  THE  REDUCED  MASS  MATRIX 


In  general,  it  may  be  expected  that  the  algebraic  sum  of  all 
the  elements  of  a  reduced  mass  matrix  from  system  identifica¬ 
tion  will  approximate  the  gross  weight  of  the  aircraft.  Due 
partly  to  restraints,  the  sum  of  the  elements  should  not 
exactly  equal  the  gross  weight,  because  masses  at  elastic 
restraints  do  not  act  as  if  they  were  ungrounded.  Masses  at 
pinned  joints  to  ground  do  not  even  figure  in  the  mass  matrix 
because  they  do  not  move. 

Individual  mass  elements  cannot  be  interpreted  as  reflecting 
lumped  physical  weights  at  their  assigned  locations.  The 
elements  of  any  reduced  mass  matrix  represent  the  inertial, 
as  opposed  to  elastic  and  damping,  dynamic  effects  of  the 
two  (for  off-diagonal)  degrees  of  freedom  with  which  they  are 
associated  in  an  actual  system  having  many  more  degrees  of 
freedom  than  the  model.  The  off-diagonal  terms  in  a  reduced 
mass  matrix  will  usually  be  large  and  sometimes  negative. 

The  matrix  will  usually  be  fully  populated. 

The  identified  mass  and  stiffness  matrices  can  be  used  to 
draw  a  dynamic  circuit  of  the  helicopter  and,  if  any  one  were 
interested,  it  would  be  possible  to  construct  an  actual  spring- 
mass  system  (utilizing  both  positive  and  negative  springs  and 
moments  of  inertia)  which  would  be  an  exact  physical  duplica¬ 
tion  of  the  identified  model,  element  by  element,  and  would 
have  the  same  natural  frequencies  and  modal  eigenvectors  as 
the  helicopter;  but  it  would  not  "look"  like  a  helicopter. 
Neither  negative  spring  rates  nor  negative  off-diagonal 
masses  are  physically  unrealizable;  the  former  are  used  by 


14 


Lockheed  in  its  control  system  and  by  all  light-switch  manu¬ 
facturers,  the  latter  are  the  essential  part  of  the  dynamic 
antiresonant  vibration  isolator. 

The  objective  is  not  to  identify  a  system  which  "looks"  like 
a  helicopter  but  one  which  "performs"  like  a  helicopter  under 
various  dynamic  loadings.  The  physical  interpretation  of  the 
ij-th  element  of  the  identified  mass  matrix,  for  example, 
is  that  the  helicopter  will  generally  exhibit  a  partial 
derivative  of  a  force  at  i  with  respect  to  a  response  at  j 
which  has  an  effective*  mass  component  that  is  the  ij-th 
element  of  the  identified  mass  matrix  (similarly  for  the 
stiffness  and  damping  matrices)  . 

It  is  immaterial  in  the  identification  whether  there  are  as 
many  points  on  the  structure  as  there  are  degrees  of  freedom 
in  the  model,  or  if  up  to  three  degrees  of  freedom  (in 
orthogonal  directions)  occur  at  any  one  point.  It  is  impor¬ 
tant  only  that  elements  in  the  motion  vector  have  the  proper¬ 
ties  of  generalized  coordinates  for  the  holonomic  model 
considered.  An  identified  reduced  model  in  which  some  of 
the  displacement  elements  represent  the  orthogonal  Cartesian 
or  polar  coordinates  of  a  given  structural  point  would  look 
much  like  an  identified  model  of  a  similar  system  with  parallel 
coordinates  of  separate  points. 

The  impedance  matrix,  of  which  the  mass  and  stiffness  matrices 
are  terms,  of  a  mathematical  model  of  a  larger  system  is  a 
function  of  the  size  of  the  model,  and  the  terms  must  reflect 
this.  It  was  found  that  frequency-independent  mass,  stiff¬ 
ness  and  damping  matrices  as  described  can  accurately  reflect 
the  responses  of  a  continuous  structure  over  a  finite  spectrum 
by  approximating  a  lambda  matrix  the  inverse  of  which  very 
closely  approximates  the  mobility.  The  spectral  mobility 
matrix,  even  of  an  order  that  equals  the  number  of  degrees 
of  freedom  in  the  structure,  cannot  be  expressed  as  a  lambda 
matrix. 


*Not  to  be  confused  with  the  formal  definition  of  "Effective 
Mass"  as 


ME 


tm]  4^ 

jki  5  *ji*ki 
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THE  REDUCED  MASS  MATRIX 


Consider  the  actual  structure  to  consist  of  an  infinite 
number  of  degrees  of  freedom  of  which  R  degrees  of  freedom 
are  retained  in  the  model.  The  mobility 
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(31) 


The  model  impedance  is  defined  as  the  inverse  of  the  mobility 
matrix  in  the  R  degrees  of  freedom: 


[Zm]  "  ^RR1  =  IZRR]  “  IZRE5  [ZW  tZER]  =  *KmJ  ”  “  lMm] 

(32) 

The  stiffness  of  the  model,  [K^]  ,  is  the  inverse  of  the  RxR 
influence  coefficients: 

*KmJ  =  ^CRR]  =  ^KRR^  “  [KRE] [KEE]  [KER^  (33) 

From  Equation  (31)  , 
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Substitute  into  Equation  (32) 
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Substitute  Equation  (33)  •  Then 
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-  <«2tKmrl[MEEj|  -1  [K^]-1!^] 


I1’ 

He1)  (l 


tMmJ  *  [fW  +  <KRE1[KEE1'1[MEEI  ([I! 

-  He’He1)  He^Hr1  (34) 

Equation  (34)  is  the  "exact"  RxR  reduced  mass  matrix  of  a 
system  with  an  infinite  number  of  degrees  of  freedom.  Note 
that  [M  ]  is  not  diagonal  and  is  a  function  of  forcing 
frequency. 


The  frequency  dependency  of  the  "exact”  reduced  mass  matrix 
simply  reflects  the  fact  that  R  linear  differential  equations 
with  constant  coefficients  cannot  contain  enough  information 
to  exactly  reflect  the  action  of  an  infinite  number  of  degrees 
of  freedom  over  a  spectrum  containing  R  modes.  The  frequency 
dependency  makes  it  impractical  to  use  this  in  a  linear 
engineering  mathematical  model. 


The  "Consistent  Mass  Matrix"  (Reference  3) ,  often  used  in 
finite-element  dynamics  work,  is  also  based  on  a  model  stiff¬ 
ness  matrix  [Km]  being  the  inverse  of  the  RxR  influence  co¬ 
efficient  matrix: 

[Kml  =  [CRR]  =  [Krr]  -  [Kj^,]  [ke£]  [rer1  • 

The  kinetic  energy  of  the  structure  is  set  equal  to  the 
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kinetic  energy  of  the  model: 


+  {yE}  [M^HyjJ  (yR>  IMmHyRl 


(35) 


It  is  implicitly  assumed,  however,  that  the  inertial  forces 
occur  only  along  the  R  generalized  coordinates ,  giving 


[Wrr1 

ikre'1 

yR  j 

Iker^ 

Ikee  \ 

lyE! 

I  0 

which  is  clearly  not  the  case  but  from  which  it  follows  that 
{yEl  =  -  [ree]  ^ [KRR] {yR}  in  sinusoidal  vibration.  Sub¬ 
stituting  the  above  in  Equation  (35)  gives 

^R^MRr]  +  (fKER^  ^KEE^  ^E^^EE^  ^KER^  ^YR^ 

=  tyR>TtMmif^R}  (36) 

and  the  "Consistent  Mass  Matrix"  is  given  by 

[Mm’  -  [MRR]  +  ‘KER'T'KEEI'T'MEE’>KEE’'ltKER]  <37> 


This  matrix  is  nondiagonal,  like  the  "exact"  reduced  mass 
matrix,  and  has  the  advantage  of  being  independent  of  fre¬ 
quency.  However,  comparison  of  Equation  (37)  with  Equation 
(34) shows  that  the  "Consistent  Mass  Matrix"  reduces  to  the 
"exact"  reduced  mass  matrix  only  at  zero  frequency;  that  is, 
in  the  static  condition.  As  the  frequency  increases,  the 
"Consistent  Mass  Matrix"  yields  increasingly  erroneous 
results . 
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The  reduced  mass  matrix  in  system  identification  is,  like 
the  others,  nondiagonal  and  related  to  a  model  stiffness 
matrix  which  is  the  inverse  of  the  RxR  influence  coefficients 
(as  represented  by  the  first  R  modes,  which  is  accurate 
beyond  direct  measurement  capability  by  many  orders  of 
magnitude) ?  but  the  system  identification  reduced  mass  matrix 
also  is  independent  of  frequency  and  is  exact  at  all  the 
natural  frequencies  of  the  model,  which  are  the  first  R 
natural  frequencies  of  the  helicopter.  The  system  iden¬ 
tification  mass  matrix  is  given  by 

[Mm]  =  IMrr]  +  t$RR]  [®ER1  £MEE]  [$ER]  [<tRR] 

or  [MJ  ^  [Crr]_1[$rr][-T]  ^RR^"1  ver^  nearly.  (39) 

qr 

At  the  r-th  natural  frequency , 

[CRR]  1  [[MRR]  +  tKRE]  [KEE]  (jjT  [I1 

r 

[CEE]  lMEE^  fKEE]  [KER]  ^RR1  [MRrJ  t<J>Rr*  =  ^Rr*  Q  2 

y r 

(40) 

exactly.  Note  in  Equation  (38)  that  the  reduced  system 
identification  mass  matrix  is  expressed  in  terms  of  the 
modal  eigenvectors  of  the  first  R  modes  only  but  includes 
all  the  masses  of  the  actual  helicopter. 

Alterations  in  masses  on  the  R  generalized  coordinates  which 
do  not  affect  the  modal  eigenvectors  are,  as  seen  from 
Equation  (38) ,  exactly  represented.  Such  alterations  can 
substantially  change  natural  frequencies  and  responses. 

Other  types  of  changes  which  do  alter  the  modal  eigenvectors 
may  or  may  not  be  accurately  reflected  in  the  model  response 
depending  on  the  degree  of  eigenvector  effects-  a  limit  which 
has  not  been  algebraically  defined  for  any  mathematical  model, 
whether  from  intuitive  analysis  or  system  identification. 


That  such  a  limit  should  somewhere  exist  is  a  practical 
engineering  fact.  One  cannot  expect  to  obtain  the  equations 
of  a  sweet  pea  on  a  rubber  band,  then  attach  it  to  the  Golden 
Gate  bridge  and  expect  to  find  the  dynamic  response  of  the 
bridge  (the  reverse,  incidentally,  is  equally  impractical). 
Prudence  marks  the  boundary  between  utility  and  uselessness. 

INFORMATION  LOSS  IN  MATRIX  INVERSION 


It  is  inevitable  that  there  will  be  a  loss  in  information 
in  numerically  obtaining  the  response  matrix  from  any  mathe¬ 
matical  model,  or  in  obtaining  the  mathematical  model  from 
responses,  even  if  no  deliberate  error  is  introduced. 

The  following  is  a  slight  modification  of  a  derivation  by 
Rosanoff  and  Ginsburg  (Reference  4).  Consider  the  equation 

[A) {x}  =  (b)  (41) 


in  which  [A]  is  a  real  symmetric  nonsingular  matrix.  Because 
we  calculate  with  numbers  which  have  a  finite  number  of 
digits,  we  actually  solve  the  equation  4 

([A]  -  [E])  (x  +  Sx)  =  (b)  (42) 

where  [E]  is  an  "error"  matrix.  Premultiplying  both  sides  of 
Equation  (42)  by  (A]-1  and  substituting  [A] “1(b)  =  {x}  gives 


([I]  -  [A]"1[E])  (x  +  <5x)  =  (x) 


(43) 


or 


{<$x}  =  p]  -  [A]-1[E])  _1  -  [I]J{x}  (44) 

Take  the  norm  (see  References  5  and  6,  for  example)  of  both 
sides : 


liuxiu  =  ||  (m  - 


But  the  norm  of  the  product  of  a  matrix  and  a  vector  is  less 
than  the  product  of  the  matrix  norm  and  the  consistent 
vector  norm: 


1 1  Mil 


or 


(6x) 


TxT 


[([I]  -  [A]"1  [E])  -1  -  [I]J  |HUx}||  (45) 

([I]  -  [A] _1  [E])  _1  -  [I]  |  | 
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EwnfnmM 


Assume  that  ||[A]“^[E]||  <  1.  From  Faddeeva  (Reference  5), 
it  is  well  known  that 


[A]  :[E])  -1-  [I]  +  ([A] _1  [E]) 


2 


<  1 | [A] _1 [E] 1 |k+1 
“  1  -  M  [A]_1[E]  |  | 

if  | | [A]  ^[E] | |  <  1.  Setting  k  =  0  gives 


(46) 


||  ([I]  -  [A] _  1  [E ])  _1  -  (I]||  <  I  I  [A]"1  [E]  |,| -  (4?) 

V  '  ~  1  -  |  |  tA]  1[E]  |  | 

or 

I  I  {<Sx>  |  |  <  I  I  [A]  [E]|  i  <  I  I  [A]  1  I  *  I  I  [E]  |  | 

|  |  (x>  |  )  “  1  -  ||  [A]_1  [E]  |  |  “  1  -  |  |  [A]  _1 1  |  *  |  |  [E]  |  | 


which  is  identical  to  the  result  obtained  by  Rosanoff  and 
Ginsburg . 

Il»x)||  <  |  |  [A]  |  | -|  I  [Al-1!  | -|  I  [E)  |  l/H  [A]  |  |  .  V 

lUxlll  -  1  -  |  |  tA]  I  I  -  I  I  tA]_1  M  -  |  I  tE]  I  I/i  I  [A]  1  -  knJt 


(48) 


where,  following  Rosanoff  et  al,  k  is  defined  as  a  con¬ 
ditioning  number  and  t  as  a  relative  error: 


kn  =  I  I  [Ml  I 'I  I  [Al"1 1  I 

l  =  | |  [E]  |  |/| |  [A)  |  | 


Taking  the  number  of  digits  in  the  arithmetic  as  the 

reciprocal  of  Equation  (48)  gives  an  estimate  of  the  number 
of  significant  digits  p. 


p  =  log 


101 


-  log 


10 


<$x  I 


log 


10 (1  "  V>  +  log 


10  l  log 


10kn 


but,  assuming  1>>  kn  i,  this  estimate  may  be  written 
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(49) 


P  "  10?10  r  -  Io9l0kn 

Thus,  as  shown  in  Reference  4,  the  number  of  information 
digits  q  lost  in  inverting  [A]  is  approximately 

q  =  log10kn  =  loglo!  IIA1IHI  [A]"1!  I  (50) 

This  is  true  for  any  norm.  However,  the  norm  of  a  symmetrical 
positive  definite  matrix,  subordinate  to  the  Euclidian  vector, 
is  the  maximum  eigenvalue ;  and  the  maximum  eigenvalue  of  the 
inverse  is  the  reciprocal  of  the  minimum  eigenvalue  of  the 
matrix.  Substituting  this  norm  of  [A]  into  Equation  (50) 
gives  the  lost  digits  estimate. 


X  (A) 

i  max 

'J  ■  Io9io  my—  (51) 

min 

To  illustrate  the  immense  practical  importance  of  this, 
consider  as  an  example  a  matrix  having 


X  (A) 

xTaT 


max 

min 


1.7k  x  103 


This  is  the  ratio  of  natural  frequencies  in  the  20x20  specimen 
of  the  helicopter  used  in  this  contract.  The  IBM  360  uses 
six  hexadecimal  places  resulting  in  166-1  or  16777215  as  the 
largest  decimal  mantissa  in  single  precision.  The  inversion 
of  the  k~-*-m  matrix  with  single  precision  on  the  computer  re¬ 
sults  in  an  inverse  having  (estimated)  loqj016777215  - 
logigl.72  x  103  =  3.99  significant  digits.  In  other  words, 
even  starting  with  eight  decimal  places  in  floating  point, 
we  end  up  with  approximately  four  decimal  places  of  informa¬ 
tion  in  the  inverse. 

It  is  absolutely  essential  when  dealing  with  test  data 
matrices  which  will  be  inverted  that  the  ratio  of  the 
extreme  eigenvalues  be  minimized.  Otherwise,  all  the 
physical  information  in  the  matrix  is  likely  to  be  destroyed 
in  the  inversion,  leaving  meaningless  numbers.  Test  data 
has  few  enough  significant  figures  of  information  to  begin 
with . 
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HOW  TO  MINIMIZE  INFORMATION  LOSS 


A  major  step  in  this  system  identification  process  is  the 
determination  of  the  { y }  and  {<j>}  vectors  by  iteration.  The 
matrix  involved  is  the  product  of  a  mobility  matrix  and  the 
inverse  of  another  mobility  matrix.  This  inverse  presents 
a  serious  danger  of  information  loss. 

To  minimize  the  extraneous  information  content  of  modes 
higher  than  the  order  of  the  matrix,  which  amounts  to  noise, 
and  to  narrow  the  spread  of  modal  mobilities,  the  matrix  to 
be  inverted  was  made  the  sum  of  the  dissipative  (e.g.,["I]) 
matrices  measured  near  each  natural  frequency. 


Each  dissipative  mobility  matrix  has  a  high  information  con¬ 
tent  about  the  dominant  mode  and  very  little  information 
about  other  modes.  This  minimizes  pollution  by  unwanted 
modes  but  results  in  a  very  poorly  conditioned  matrix.  For 
example,  the  10x10  imaginary  acceleration  mobility  of  a 
typical  helicopter  measured  at  3  Hz  has  an  extreme  modal 
mobility  ratio  of  106.  However,  the  sum  of  mobility  matrices 
over  the  frequency  range  is  a  matrix  having  the  same  modal 
vectors  as  a  mobility  matrix  at  any  one  frequency. 


] 

P 


N 

E 

i-1 


“*i 

Y  * 
to  i 

P 


41>4i)T 


(52) 


E  tij] 

(0 


n 

E 

i=l 


5  i*?U1H*1>T 

to 


(53) 


Therefore,  Equation  (53)  can  be  used  as  one  of  the  matrices 
in  the  modal  eigenvector  equations 


=  A{y}i  (54) 


[Y1]  [EY*I]“1{(J)}. 

0)  to  l 

to 


ot  { 4> }  i 


(55) 


The  range  of  values  from  the  maximum  to  the  minimum  in 
•**I  .  ”*T 

EY  ^  is  very  small  compared  to  the  range  of  Y^. 

to 
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If  E [Y*]  or  E [YR]  is  used  in  place  of  ElY1]  it  is  necessary 
0)  0)  0)  w 
to  normalize  each  of  the  matrices  in  the  sum  because  the 
displacement  and  velocity  mobilities  decrease  in  magnitude 
with  increased  frequency.  Normalization  on  the  root  mean 
square  of  the  matrix  elements  and  on  the  largest  element 
absolute  value  were  both  investigated  experimentally. 
Normalization  on  the  RMS  gave  results  about  as  satisfactory 
as  those  from  acceleration  mobility  and  is  preferred  over 
normalization  on  the  largest  element,  as  the  latter  is  sen¬ 
sitive  to  errors  in  one  term  which  could  throw  off  the  entire 
matrix.  However,  the  differences  in  results,  while  evident, 
were  not  dramatic. 

The  log^Q  of  the  ratio  of  the  maximum  gY  ^  to  the  minimum 

was  generally  about  .75  for  the  5x5  models  in  these  experi¬ 
ments  and  generally  around  1.8  for  the  15  x  15  models.  The 
5x5  models  performed  excellently  but  the  15  x  15  models 
performed  -apriciously . 

"*I 

If  the  engineer  could  normalize  so  that  the  matrix  [EY 
is  unity,  information  would  still  be  lost  in  the  w 
inversion  but  certainly  less  information  than  if.-*j 
the  ratio  of  extreme  values  is  very  high.  The  EY  ^  terms 

"I 

are  not  the  eigenvalues  of  E •  The  only  matrix  which  has 

a  unit  eigenvalue  matrix  is  the  unit  matrix  itself;  it  follows 
therefore  that  some  information  is  always  lost  in  the  numeri¬ 
cal  inversion  of  any  matrix  other  than  unity. 

The  matrix  we  wish  to  invert  is 

£[**]  =  !*i  (5«) 

Express  the  modal  vector  matrix  in  terms  of  its  own  eigen¬ 
vectors  [J]  and  its  own  eigenvalues  (that  is,  [J]  is  the 
eigenvector  matrix  of  the  eigenvector  matrix  of  [tk ]  “ 1  [m J)  . 

[*]  =  IJ]  [J]"1  (57) 

Substitute  Equation  (57)  into  Equation  (56). 
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meuu 


E[Y *]  =  [J]  tx  J  [J]-TfX  g  [J]T 

U)  v  v 


invert , 


"I  -1  -T 

Z  [Y1]  1  =  [J]  1 


10 


<o 


LX*J 


'J,T  fej  IJ,k]IJ,‘l 


(58) 


(59) 


The  only  operation  on  the  eigenvectors  [J]  between  Equation 
(58)  and  Equation  (59)  was  to  change  relative  positions;  all 
the  inversions  were  of  diagonal  matrices.  As  a  diagonal 
matrix  is  a  matrix  of  its  own  eigenvalues,  having  the  unit 
matrix  for  eigenvectors ,  the  central  term  may  be  treated 
rigorously  as  an  eigenvalue  matrix.  The  matrix  of  Equation 
(58)  may  be  substituted  for  [A]  in  Equation  (41) ,  and  in 
Equation  (50),  we  can  consider [a] as  the  product  of  the  three 
matrices  of  Equation  (56). 


E  [Y1]  =[A]=  [$]^Y*Jj[$]T  (60) 

<0  col 

CO 

It  is  well  known  that 

lllAlll  i  llt*]IMIP:Y**J[4]T||  (61) 

and  that 

MfEY*^[4]T||  £  MttY*iIj|MlmTll  (62) 

Therefore 

1 1  [a] 1 1  <  MmlHI  tfY*^j|MlmTll  (63) 

At  this  point  we  wish  to  substitute  eigenvalues,  but  [4>]  is 
not  symmetric  so  |  |  [4>3  |  |  ; /  |max  X^l  .  Rather,  |  |  [4>]  |  |  £ 

T 

|max  X^|.  Consider,  therefore,  the  eigenvalues  Xb  of  [4>]  [4>] 
which  is  symmetrical. 

[$]Tm  =  [L]fXb^[L]_1  =  [L]fXbg(L]T  (64) 

T 

where  [L]  is  the  orthogonal  matrix  of  eigenvectors  of  [4>]  [41]  . 


I  I  I  =  |  max/X^| 


(65) 


f 


j 
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X. 


Substitute  Equation  (65)  into  Equation  (63)  . 


[A]  |  |  <  [max  Ab| • (max  EY*J | 


'66) 


Using  Equation  (51)  ,  the  number  of  digits  lost  in  inverting 

E [Y  ]  is  approximated  by 
to  a 

ImaxA,  I  •  Imax  EY  .  I 

_  ~  1  b1  1  oil1 

q  ~  log,  h  - vr*T~ 

|minA.  |  •  |min  EY  .  | 

1  b 1  1  0)1 1 


=  log10 


ImaxA,  | 

- —  + 

|minAb| 


|  max  EY  ? j 
1  0)1  1 

"  "*T~ 

mm  EY  . 


(67) 


r-AgJ  would  equal  a  scalar  times  the  unit  matrix  only  if  the 

modal  vectors  {<}>}  were  orthogonal  (i.e.,  {c^}  { . }  =  o)  ,  a 

condition  which  could  occur  only  in  the  academic  cases  of 
uniform  mass:  [m]  =  [m][I].  in  this  case,  the  loss  of 
information  digits  would  be  indicated  by 

;*i, 


max  Y 


(01 


q  ~  logio  -  T  ;»i 

lmln  Y„n 


(68) 


0)1 


and  only  in  this  case  could  zero  information  loss  be  achieved 
by  normalizing  the  matrices  such  that  |max  y*?|/|min  y*? |  =  1 

But  the  case  is  trivial,  for  if  it  were  true,  an  inversion 
would  be  unnecessary  as  $  would  be  the  eigenvector  matrix  of 

E  [YJ  . 


0) 


co¬ 


if  the  mass  distribution  is  not  uniform  diagonal  but  the 
engineer  could  so  normalize  the  matrices  in  the  summation 

50  that  Imax  l/lmin  ??|  -  1,  14  is  seen  £rom  Ration  (67) 
1  0)1  1  1  0)1 

that  there  would  still  be  a  loss  of  information  digits 
approximated  by 
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q 


[max  A  | 

log,  A  - 

°  |n»in  Ab| 


(69) 


The  ratio  jmax  Ab|/|min  A^l  increases  with  the  order  of  the 

mobility  matrix;  that  is,  with  the  number  of  degrees  of 
freedom  of  the  model.  It  follows,  therefore,  that  there  is 
an  upper  limit  to  the  size  of  a  physically  meaningful  reduced 
complete  model  regardless  of  normalization  of  the  matrices  in 
the  summation. 


As  a  crude  "rule  of  thumb",  Figure  1  shows  the  trend  in 


NUMBER  OF  MODEL  DEGREES  OF  FREEDOM 
•  •  .  •  *"T 

Figure  1.  Reliability  of  the  Inversion  of  £  [y  ]  . 

u  w 

It  is  seen  that  the  reliability  of  the  calculation  becomes 
questionable  above  10  or  15  degrees  of  freedom.  This  does  not 
mean  that  accurate  identifications  cannot  be  made  using  the 
iterative  step  for  modes  of,  say,  20  degrees  of  freedom  but, 
rather,  that  any  one  calculation  has  a  higher  probability  of 
failure . 

In  passing,  it  should  be  noted  that  the  treatment  of  bounds 
using  matrix  norms,  as  above,  opens  up  some  highly  promising 
avenues  of  research  on  the  reliability  of  many  helicopter 


NX 
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theoretical  calculations  as  well  as  on  the  reliability  of  the 
processing  of  test  data  in  general.  Whether,  for  example, 
some  of  the  conventional  methods  of  processing  strain  gage 
data  yield  physically  meaningful  results  is  open  to  question 
in  the  light  of  the  above  method. 

WHEN  A  CALCULATION  FAILS 

*’T  **  Y  —  1 

The  most  common  mode  of  failure  of  iteration  on  [Y  ]  (E[Y  ]) 

0) 

is  catastrophic,  producing  such  absurd  values  for  one  or  more 
generalized  masses  as  negative  numbers  or  unusually  large 
numbers.  This  is  signified  also  by  a  very  large  number  of 
iterations  required  for  convergence  on  one  or  more  modes . 
Failures  almost  never  occur  with  small  number  of  degrees  of 
freedom  (e.g.,  five), and  an  identification  which  is  quite 
accurate  with  one  seed  may,  in  the  larger  models,  diverge  with 
another  seed. 

This  phenomenon  results  from  the  fact  that  the  significant 
effect  of  error  is  not  insidious  accumulation  of  inaccuracies 
in  the  generalized  masses  but,  rather,  information  destruction 
in  the  inversion.  Fortunately,  it  is  usually  very  obvious  to 
the  engineer  when  an  identification  fails  on  the  computer,  and 
corrective  measures  may  be  taken  without  rerunning  the  test 
on  the  helicopter. 

A  most  obvious  and  effective  corrective  measure  is  to  elimi¬ 
nate  one  or  more  of  the  degrees  of  freedom.  This  can  be  done 
on  the  computer  ,  as  the  program  is  written  so  that  the  system 
may  be  instructed  to  select  any  of  the  available  data  which 
is  in  digital  form  on  tape.  The  size  of  the  model  and  the 
number  of  modes  covered  are  consequently  reduced.  It  is 
possible  also  to  eliminate  any  mode,  not  just  the  highest, 
if  it  appears  that  a  certain  mode  contributes  little  in¬ 
formation  -  a  local  resonance,  for  example,  in  which  only  a 
small  portion  of  the  helicopter  is  significantly  responding. 

The  computer  experiments  included  a  local  resonance  in  the 
form  of  a  mode  in  which  only  the  most  forward  station  showed 
substantial  motion.  When  this  station  was  not  included 
in  the  identification,  but  the  local  resonance  associated 
with  its  movement  was  included,  then,  as  expected,  there 
were  evidences  of  failure  in  generalized  mass  calculation. 

The  computer  was  attempting  to  identify  a  natural  mode  for 
which  the  input  mobility  data  showed  a  largely  nonresponding 
helicopter.  This  situation  would  be  detected  from  the 
mobility  plots  before  committing  the  data,  as  it  is 
very  apparent  in  the  dissipative  mobility  spectra.  The 
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ability  to  handle  local  resonances,  or  dispose  of  them  when 
required,  is  important  to  a  practical  identification  because 
all  real  structures  have  them.  In  fact,  as  the  number  of 
degrees  of  freedom  of  a  simulated  structure  are  realistically 
increased,  the  modal  density  usually  increases  more  rapidly 
than  the  simple  mathematics  of  uniform  chains  would  lead  one 
to  believe.  When  that  degree  of  freedom  which  is  the  pre¬ 
dominant  motion  of  a  local  resonance  is  eliminated,  the  mode 
it  causes  should  be  eliminated  also,*  the  mobility  spectra 
plots  for  the  included  degrees  of  freedom  would  indicate  this 
by  an  insignificant  peak. 

THE  r  MATRIX  AND  MODAL  PARAMETERS 


The  dominant  modal  vector  at  frequency  wi,  near  the  i-th 
natural  frequency  ^i,  is  given  by  Equation  (55)  and  the  i-th 
gamma  vector  by  iteration  on  the  transpose,  Equation  (54)  . 
The  modal  mobility  is  obtained  from 


{yJT  [Y  .  ]{y.  } 

(ITR)  wi  1 


Y*1 

(ITR)  la)i 


(70) 


where  (y.)  is  the  vector  from  iteration  (Equation  54). 

1  (ITR) 

It  is  impractical  to  attempt  the  calculation  using  {y.}  from 
because  of  information  loss  in  the  inversion,  as  shown 
in  Equation  (59)  .  The  dominant  mode  is  the  only  one  used,  of 
course,  as  there  is  negligible  information  content  in  ^"1  ^ 
about  modes  other  than  the  i-th.  Therefore, 


[Y«i> £  yii{*i,{*i,T  s  YIL{*i>  {+<>' 


(ITR) 


(ITR) 


(71) 


and  (y.}T  {<(>.  }  =  1  is  forced. 

(ITR)  1  (ITR) 

A  peculiar  situation  often  occurred  when  a  calculation 
diverged:  it  was  noticed  that  the  natural  frequency  of  the 
"bad"  mode  was  usually  identified  with  great  accuracy  although 
the  calculated  generalized  mass  was  absurd,  often  negative, 
and  negative  calculated  values  of  ^*R  often  occurred. 


X 


X 

X 
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*  *R 

The  key  here  is  the  occurrence  of  negative  values  cf  Y 

•  R  1Wi 

Ideally,  [y  ]  is  a  positive  definite  matrix  and  cannot, 

theoretically,  be  negative  definite  on  grounds  that  it 

represents  the  dissipation,  not  a  source,  of  energy.  For 

any  positive  definite  matrix  B,  {x)T[B]{x}  is  a  positive 

number  regardless  of  the  choice  of  the  vector  (x).  The  fault 

for  negative  values  of  **R  ,  which  are  physically  impossible, 

*ia). 

1 

cannot  therefore  be  laid  solely  to  {y},  and  therefore  to  the 
loss  of  information  in  the  inverse  of 


<V 

(iRl 


T  ,:.R, 


£ 

a) 


r*I,  because 
1  uJ  ' 


[Yx]{y^}  must  be  positive  even  for  arbitrary  {y^}  if 
is,  as  it  is  supposed  to  be,  positive  definite.  We 


forced  to  conclude  that  numerical  errors  can  act  in  such 

i.Rj  not  positive  definite. 


are 

a 


way  as  to  make 

The  mobility  [Y*'  ]  is  very  nearly  equal  to  the  positive  semi- 


[r 

■R 


definite  matrix 
positive.  Then 


0). 


.  •  *R 

in  which  Y. 

10)  . 


is  necessarily 


(72) 


But  {y.}  and  {<(>.  }  are  composed  of  real  numbers,  as  opposed  to 

^  ^  Ijl 

imaginary  or  complex  numbers,  which  makes  {y.  >  {<}>.}  real  and 

( (y.  n<|>  • )) z  real  and  positive  even  for  arbitrary  elements  in 
^  ^  •  R 

{y}.  The  dominance  of  [Y  .]  by  one  mode  is  therefore  not  a 

.  **R 

cause  of  calculating  negative  values  of  Y^u  . 

i 

•  *R 

The  calculation  of  absurd  values  of  Y.  is  nevertheless  due 

10)  ••  j 

mainly  to  information  loss  in  inverting  EtY^l  (or  other 

normalized  mobility  matrices  having  similar  properties),  which 
results  in  poor  eigenvectors  in  the  iteration.  Examination 
of  the  computer  experiments  shows  that  errors  in  the  [YR]  or 
[Y1]  matrices  are  not  sufficient  to  cause  as  erratic  results 
as  have  sometimes  been  observed  if  the  { y }  vectors  in 
{y}T[Y]{y}  are  accurate.  In  the  "bad"  cases,  the  {y}  vectors 
from  iteration  are  invariably  very  bad.  The  reason  for  the 
occasional  negative  calculated  values  of  **R  is,  in  part, 
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that  errors  in  [Y]  can  cause  the  matrix  to  not  be  positive 
definite.  For  example,  in  Computer  Experiment  188  a  nine- 
point  identification  with  error  yielded  good  results  but  the 
same  identification  with  a  different  seed  (Computer  Experi¬ 
ment  184)  gave  poor  results  which  included  negative  **R  for 

Yito 

the  seventh  mode.  The  errors  by  chance  happened  to  act  in 
such  a  way  in  Experiment  184  that  excessive  information  was 
lost  in  the  inverse,  as  indicated  by  iterations  that  failed 
to  converge.  The  principal  minor  associated  with  the  eighth 
and  ninth  positions  in  mobilities  dominated  by  the  seventh 
mode  was  found  to  be  negative  in  the  bad  case  (Experiment  184), 
due  to  a  peculiar  accumulation  of  random  errors,  which,  of 
course,  meant  that  the  mobility  was  no  longer  positive  def¬ 
inite,  as  in  pure  theory,  and  could  give  negative  values  of 
T 

{y}.  [Y  Hy  •  } .  However,  precise  {y}  vectors  would  not  have 
i  a) .  l 

X  •  *R 

caused  the  negative  values  of  Y.  even  with  [Y]  not  being 
positive  semi definite. 

*R 

Calculation  of  physically  meaningless  values  of  Y.  ,  and 

therefore  of  M*,  is  caused  primarily  by  information  loss 
in  inversion. 

The  reason  for  fairly  accurate  identifications  of  natural 
.frequencies  even  when  the  generalized  mass  identifications 
are  poor  lies  in  the  fact  that  co.  and  cok  in  Equation  (22) 
are  taken  near  fh  ;  therefore,  ^ 


~  1  (73) 


and 


2 

oo.oo,  =  ft.  +  .  ( 6  uj,  -  6  u) . )  -  6co.5io 

]  k  i  i  k  j  Dk 


2  ~ 

But  6wi,  <<1>>  6  oj .  so  fl.  ~  w.ojv 
K  D  1  3  * 


•  *1 

<o.Z. 

3  iio 

“i^k  — r*~ 

D  z. 

k  loo. 


•  *1 

CO,  Z. 
k  i(0j 

•  *1 

co  .Z . 

3  ico. 


co.co,  =  (f2.  -  6 co . )  ( Q .  +  6 co,  ) 
1  k  i  T  i  k 


. * j  . *1 

co.Z.  -  co,  Z. 
j  icok  k  ico. 


Ti  • 

0).  Z  .  -  (A)  .  Z  .  ‘ 

k  i(ok  3  i(0j 
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IDENTIFIED  GENERALIZED  MASSES 


Typical  generalized  mass  identifications  are  shown  in  Tables 
I  through  VI.  Note  in  Tables  I,  III  and  V  that  the  generalized 
mass  of  the  first  mode  identified  for  a  reduced  model  with  no 
experimental  error  has  always  been  less  than  the  first  mode 
generalized  mass  calculated  from  the  modal  vector  and  mass 
matrix  of  the  specimen.  This  is  not  true  of  other  inodes. 

Tables  I  and  II  show  results  of  two  different  five-point 
models.  No  outstanding  differences  between  the  models  is 
evident.  Model  9A  produced  acceptable  results,  as  shown  in 
Table  III,  for  different  distribution  of  random  error  but 
Model  9B,  as  shown  in  Table  IV,  worked  with  some  seeds  and 
failed  with  other  seeds.  The  failed  experiments  of  Table  IV, 
Computer  Experiments  168  and  184,  yielded  drastically  un¬ 
realistic  values  of  generalized  mass  for  most  of  the  modes. 

Table  V  shows  a  twelve-point  model  identification  which 
failed  only  in  the  eighth  mode.  Computer  Experiment  178  is 
identical  to  Computer  Experiment  169  except  that  in  the 
former,  the  computer  was  instructed  to  skip  the  eighth  mode 
and,  instead,  operate  on  tape  data  for  the  thirteenth  mode 
which  resulted  in  satisfactory  identification. 

Using  different  stations  for  a  twelve-point  model,  as  shown 
in  Table  VI,  produced  proper  identification  of  all  models, 
including  the  eighth,  with  various  error  distributions. 

Information  loss  in  the  inversion  of  mobility  matrices  is 
the  primary  cause  of  such  failures,  as  shown  in  Computer 
Experiments  168,  184  and  169.  The  averaging  of  mobility 
test  data,  properly  done,  would  greatly  minimize  the  chances 
of  such  identification  failures.  Test  data  averaging  is  the 
customary  practice.  These  computer  experiments  did  not  take 
advantage  of  averaging  experiments. 
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TABLE  I.  IDENTIFICATION  OF  GENERALIZED  MASSES, 
5X5  MODEL*  OF  20  X  20  SPECIMEN 


Computer 

Experiment 

Number 

152 

151 

157 

160 

182 

Random  Amp  Error 

0 

+  5% 

+5% 

+  5% 

+  5% 

0 

Bias  Amp  Error 

0 

+  5% 

+  5% 

+5% 

+  5% 

0 

Random  Phase  Error 

0 

+  1 

+3.° 

+1° 

+1° 

0 

Seed 

- 

246 

221 

195 

327 

- 

Stations 
(In.)  Mode 

Generalized  Masses 
(Lb-Sec2/In.) 

12 

22 

34 

46 


*  Model  5 A 

**  From  20  x  20  Model 


7.9910  7.3594  7.3834  7.8421  8.2330  8.5341 

4.6248  3.8247  4.5951  4.1440  4.0594  4.4491 

.4951  .4618  .4771  .4653  .4729  .4951 

1.0897  1.0372  1.0657  1.0366  1.0440  1.0872 

.6463  .5869  .6247  .6691  .6131  .6302 


TABLE 

II 

.  IDENTIFICATION  OF 
5X5  MODEL*  OF  20 

GENERALIZED  MASSES, 

X  20  SPECIMEN 

ii  1 

Computer 

Experiment 

Number 

i 

159 

170 

183 

1** 

Random  Amp  Error 

+  5% 

+  5% 

+  5% 

0 

Bias  Amp  Error 

+  5% 

+  5% 

+  5% 

0 

Random  Phase 

Error 

+  1° 

+  1° 

+1° 

- 

Seed 

221 

246 

128 

— 

Stations 

(In.) 

Mode 

Generalized 

(Lb-Sec2/I 

Masses 
n. ) 

0 

1 

7.4385 

7.3210  7 

.8179 

8.5341 

100 

2 

4.4545 

4.1824  4 

.3797 

4.4491 

200 

3 

.4724 

.4620 

.4596 

.4951 

320 

4 

1.0769 

1.0277  1 

.0233 

1.0872 

460 

5 

.6912 

.5945 

.6360 

.6302 

*  Model  5B 

**  From  20  x 

20 

Model 

.  —   . .  1 

34 
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TABLE  III.  IDENTIFICATION  OF  GENERALIZED  MASSES, 
9X9  MODEL*  OF  20  X  20  SPECIMEN 


Computer 

Experiment 


Number 

180 

156 

162 

179 

187 

1** 

Random  Amp  Error 

0 

+  5% 

+  5% 

+  5% 

+  5% 

0 

Bias  Amp 

Error 

0 

+  5% 

+  5% 

+  5% 

+  5% 

0 

Random  Phase  Error  0 

+1° 

+  1° 

l3-0 

+  1° 

0 

Seed 

- 

287 

50 

315 

492 

- 

Stations 
(In. ) 

Mode 

Generalized  Masses 
(Lb-Sec2/In.) 

0 

1 

7.9538 

7.3776 

7.5946 

8.2378 

7.4531 

8.5342 

30 

2 

4.5889 

4.1130 

4.2450 

4.6233 

4.2020 

4.4491 

100 

3 

.4938 

.4671 

.4821 

.4614 

.4656 

.4951 

160 

4 

1.0863 

1.0507 

1.0368 

1.0129 

1.0785 

1.0872 

220 

5 

.6350 

.6164 

.0044 

.6102 

.5971 

.6302 

280 

6 

.7457 

.7049 

.6983 

.7227 

.7239 

.7429 

340 

7 

1.1746 

1.1204 

1.1332 

1.1064 

1.0968 

1.1769 

400 

8 

1.5002 

1.3770 

1.4070 

1.4193 

1.4783 

1.4683 

460 

9 

.6593 

.6576 

.5507 

.6235 

.5737 

.7866 

*  Model  9 A 

**  From  20  x  20  Model 
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TABLE  IV.  IDENTIFICATION  OF  GENERALIZED  MASSES, 
9X9  MODEL*  OF  20  X  20  SPECIMEN 


Computer 

Experiment 

Number 

±tl 

188 

168  184 

1** 

Random  Amp  Error 

±5% 

+  5% 

+5%  ±5% 

0 

Bias  Amp  Error 

+  5% 

■>•5% 

+5% 

0 

Random  Phase  Error 

±1 

±1 

+  1  ±1 

- 

Seed 

287 

206 

395  619 

- 

Stations 
(In. ) 

Mode 

Generalized  Masses 
(Lb-Sec 2/In. ) 

0 

1 

7.4445 

7.5891 

7.6797  7.0969 

8.5341 

60 

2 

4.2851 

4.4084 

23.2234  4.5730 

4.5615 

120 

3 

.4741 

.4545 

.6876  .4314 

.4951 

180 

4 

1.0194 

1.0226 

28.5896  1.0968 

1.0872 

240 

5 

.6343 

.6740 

.5667  -7.9847 

.6302 

280 

6 

.7020 

.6987 

-8.5143  .5237 

.7429 

320 

7 

1.1877 

1.0711 

-.0080  .0125 

1.1769 

400 

8 

1.2510 

1.7815 

.1256  -.2199 

1.4683 

460 

9 

.9347 

.9398 

-.0159  -.0810 

.9836 

*  Model  9B 

**  From  20  x  20  Model 
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TABLE 

V. 

IDENTIFICATION 
12  X  12  MODEL* 

OF  GENERALIZED  MASSES, 

OF  20  X  20  SPECIMEN 

Computer 

Experiment 

Number 

169 

178 

1** 

Random  Amp  Error 

+  5% 

+  5% 

0 

Bias  Amp  Error 

+  5% 

+5% 

0 

Random  Phase 

Error 

+1° 

+1° 

0 

Seed 

492 

492 

_ 

Stations 

(in.) 

Mode 

Generalized  Masses 
(Lb-Sec2/In.) 

0 

1 

7.4551 

7.4629 

8.5341 

60 

2 

4.1298 

3.9789 

4.4491 

100 

3 

.4587 

.4657 

.4951 

120 

4 

1.0446 

1.0376 

1.0872 

160 

5 

.5950 

.5802 

.6302 

200 

6 

.6869 

.6975 

.7429 

240 

7 

1.2036 

1.2044 

1.2569 

280 

8 

-7.9616 

2.0521 

320 

9 

.9410 

.9118 

.9836 

370 

10 

.0425 

.0428 

.0432 

430 

11 

.1718 

.1752 

.1723 

460 

12 

1.0012 

1.0037 

1.0480 

13 

- 

.7924 

.5724 

*  Model  12D 


**  From  20  x  20  Model 


TABLE  VI. 

IDENTIFICATION  OF 
12  X  12  MODEL*  OF 

GENERALIZED  MASSES 
20  X  20  SPECIMEN 

9 

Computer 

Experiment 

Number 

150 

149 

155 

163 

1** 

Random  Amp 

Error  0 

+  5% 

+  5% 

+  5% 

0 

Bias  Amp 

Error 

0 

+  5% 

+  5% 

+5% 

0 

Random  Phase  Error  0 

+  1° 

l10 

+1° 

0 

Seed 

- 

23 

492 

87 

- 

Stations 

(In.) 

Mode 

Generalized  Masses 
(Lb-Sec2/In.) 

0 

1 

7.9718 

7.7160 

7.2917 

7.4071 

8.5342 

30 

2 

4.6071 

4.5010 

4.2722 

4.3406 

4.4491 

60 

3 

.4941 

.4640 

.4682 

.4611 

.4951 

100 

4 

1.0857 

1.0499 

1.0625 

1.0425 

1.0872 

140 

5 

.6348 

.6094 

.5958 

.5936 

.6302 

180 

6 

.7441 

.7155 

.6930 

.7097 

.7429 

220 

7 

1.1765 

1.1433 

1.1101 

1.1278 

1.1769 

260 

8 

1.4158 

1.3467 

1.3225 

1.3454 

1.4115 

300 

9 

.7808 

.7329 

.7395 

.7362 

.7866 

340 

10 

.0430 

.0419 

.0422 

.0422 

.0432 

400 

11 

,1705 

.1596 

.1665 

.1689 

.1723 

460 

12 

.9112 

.5712 

.8417 

1.0946 

1.3235 

*  Model 

12B 

**  From  20 

x  20 

Model 

RESPONSE  FROM  IDENTIFIED  MODEL 


Figures  2  through  7  portray  typical  acceleration  response 
obtained  from  the  various  models  investigated  in  the  present 
study-  In  each  instance,  the  exact  curve  was  obtained  from 
the  twenty-point  structure  with  zero  error.  Figure  2 
indicates  the  effect  of  random  number  seed  for  a  typical 
five-point  model.  Figure  3  presents  the  results  obtained 
for  one  of  the  nine-point  models  considered  in  the  investi¬ 
gation.  Figure  4  portrays  the  effect  of  random  number 
seed  on  the  twelve-point  model.  All  the  computer  experiments 
which  considered  error  used  a  +5%  random,  5%  bias  and  a  1° 
phase  error. 

Figure  5  presents  the  effect  of  model  variation  on  the 
acceleration  response.  The  models  varied  in  that  different 
spanwise  masses  were  considered.  Model  5A  utilized  stations 
0,  120,  220,  340  and  460  (inches)  whereas  model  5B  consisted 
of  stations  0,  100,  200,  320,  anu  460  (inches).  Figure  6 
presents  the  effect  of  model  for  the  nine-point  model.  The 
model  9A  consisted  of  stations  0,  30,  100,  160,  220,  2R0,  340, 
400  and  460  (inches).  Model  9B  included  stations  0,  60,  120, 
180,  240,  280,  320,  400  and  460  (inches).  The  twelve-point 
model  12B  used  stations  0,  30,  60,  100,  140,  180,  220,  260, 
300,  340,  400  and  460  (inches)  whereas  model  12E  utilized 
stations  0,  30,  60,  100,  120,  160,  200,  260,  280,  340,  400, 

460  (inches)  .  For  each  model  the  computer  experiments  were 
executed  using  the  same  random  number  seed  and  the  afore¬ 
mentioned  errors  were  incorporated. 


IDENTIFIED 
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2  10  100 

FORCING  FREQUENCY  (CPS) 

Figure  6.  Nine-Point  Model  Response , Effect  of  Model; 
Driving  Point  at  Station  1. 
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CONCLUSIONS 


1.  The  equations  of  motion  for  a  structure  may  be  obtained 
using  only  impedance- type  test  data  without  the  use  of 
an  intuitive  mathematical  model. 

2.  The  method  also  yields  the  eigenvector  or  mode  shape 
and  generalized  mass  corresponding  to  each  natural 
frequency. 

3.  The  accuracy  of  the  dynamic  response  of  a  structure 
using  impedance- type  experimental  data  is  not  dependent 
on  the  accuracy  of  the  test  measurements,  provided  the 
data  is  within  the  state  of  the  measurement  art. 

4.  The  mass  matrix  assumed  for  an  intuitive  mathematical 
model  should  be  fully  populated  to  yield  accurate 
dynamic  response  results. 

5.  To  insure  minimum  information  .loss  in  the  inversion 
of  mobility  matrices,  the  avenging  of  mobility  test 
data  should  be  used  in  practice 

6.  There  is  an  upper  limit  to  the  size  of  a  physically 
meaningful  reduced  complete  model  yielding  minimum 
loss  of  information  digits.  The  present  renort 
indicates  the  maximum  to  be  a  model  of  approximately 
15  degrees  of  freedom. 
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APPENDIX 

COMPUTER  PROGRAM  DESCRIPTION 
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LIST  OF  FORTRAN  SUBROUTINES 


I 


AMP 

CINV 

ERR 

GEN 

INVRS 

ITER 

MITER 

MMPX 

Mgfe 

MCftjT 

RANDU 

RED 

Y0UT 

SYM 

MpUT2 

MMPY 

SITER 

MATAMP 

YRNRM 


Converts  mobility  from  velocity  units  to 
acceleration  as  amplitude  (in  g's)  and  phase 
angle  (in  degrees) 

Complex  inverse  of  complex  matrix 

Incorporates  measurement  errors  into 
simulated  measurements 

Generalized  function  of  form  f  Af  where  f  is 
a  vector  and  A  is  a  square  matrix 

Inverse  of  a  matrix 

Matrix  iteration  for  eigenvalues  and 
eigenvectors 

More  general  iteration  on  product  of  two 
matrices;  used  for  gamma  iteration 

Matrix  multiplication 

Calculates  complex  impedance  and  mobility 
Special  output  for  square  matrix 
Random  number  generator 
Removes  rows  and  columns  from  matrix 
Special  matrix  output 

Forms  symmetric  matrix  from  lower  triangle 
Special  output  for  nonsquare  matrix 
Matrix  multiplication 

Matrix  iteration  for  eigenvalues  and  eigenvectors 

Converts  velocity  mobility  to  amplitude  (g's) 
and  phase  (degrees) 

Performs  normalization  of  mobility  matrix  on 
absolute  value  of  largest  element  of  mobility 
matrix 
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YRRMS  Performs  normalization  of  mobility  matrix  on  root 
mean  square  value  of  mobility  matrix 

M0B2  Calculates  complex  impedance  and  mobility 

YRFREQ  Multiplies  each  velocity  mobility  matrix  by 

its  respective  frequency  to  give  acceleration 
mobility 

REIDN  Identification  of  mass  stiffness  and  damping 
matrices 
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INtXALT  SIMULATED  TEST  DATA 
AACT  20  POINT  UH2  8/19/70 

INXACT  9  POINT  NOD  EL  20  POINT  STRUCTURE  1-/11/70 
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IDENTIFIED  DAMPING  MATRIX 
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DRIVING  POINT  RESPONSE,  IMP  IN  G*  S  AND  PHASE  IN  DEGREES 
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